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Lattice Monte Carlo simulations are performed for the SU(2) Yang Mills gauge theory in the presence of 
an Abelian background with external sources to obtain information on the effective potential. The goal is to 
investigate the lowest Landau mode that, in the continuum one- loop effective potential, is the crucial mode for 
instability. It is shown that also in the lattice formulation this lowest Landau mode plays a very peculiar role, 
and it is important for the understanding of the vacuum properties. 



> 

(N 
O 
O 

iu: 



X 



1. INTRODUCTION 

To understand the IR behavior of non- Abelian 
gauge theory a non-perturbative framework is 
necessary. Therefore, lattice formulation is par- 
ticularly useful. However, the comparison be- 
tween one-loop perturbative expansion and lat- 
tice regularization might give important informa- 
tion. From this comparison we can learn about 
the trustworthiness of perturbative expansion. 
Furthermore, the loop expansion can provide in- 
dicative information for the lattice quantities. 

Despite the importance of Yang-Mills theories, 
the complete solution of non- Abelian gauge the- 
ories has yet to be found. In order to gain a bet- 
ter understanding of these theories, a necessary 
first step is the study of their vacuum structures. 
Nevertheless, the vacuum structure of such theo- 
ries has yet to be understood even in the simplest 
case, i.e. SU(2) without matter. Moreover, the 
infrared properties of such theories must be stud- 
ied systematically if we want to have some clue 
on confinement. 

Many authors have studied one-loop effective 
potential for SU(2) and the possible consequences 
jl||^ . In addition, several Monte Carlo simula- 
tions have already been done with a wild range 
of techniques in 3 and 4 Dimensions . 
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2. ONE-LOOP APPROACH 

A powerful method to investigate the proper- 
ties of Yang-Mills theories is to compute of the 
effective potential in the background gauge. This 
method has been discussed extensively in the lit- 
erature [Q. The technique is to split the gauge 
field into a background and a quantized field 
as 



A' 



A. 



(1) 



and subsequently performing a loop expansion. 
This manifestly gauge invariant scheme is based 
on the observation that the loop expansion corre- 
sponds to an expansion in the parameter h which 
multiplies the entire action. Hence, a shift of the 
fields or a redefinition of the division of the La- 
grangian into free and interacting parts can be 
performed at any finite order of the loop expan- 
sion without violating the gauge invariance. Nev- 
ertheless, there are several subtleties on the ex- 
ponentiation of the gauge constraint and in the 
ghosts contributions to the finite part of the ef- 
fective action 

We want to investigate the situation where A 
generates a static constant chromomagnetic field. 
A possible choice for A^^ is the so called Abelian 
background: 



Al = \H5''^{x5^2-y5^i) . 



(2) 
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With this choice, the well known Savvidy result 
for the one-loop effective potential is obtained 

The remarkable feature of this expression is that 
it exhibits a minimum for H different from zero. 
However, as Nielsen and Olesen realized, due 
to the imaginary part of the effective potential, 
this minimum has unstable modes. Since the 
above mentioned preliminary studies were done, 
the property of this non-trivial vacuum has been 
thoroughly investigated. In particular a scenario, 
the so-called "Copenhagen vacuum" was pro- 
posed. 

However, in a strong field configuration a 
perturbative analysis is unreliable, and unsta- 
ble configurations can only be analyzed by non- 
perturbative methods. Therefore, the only pos- 
sible technique presently available to tackle this 
problem is the one based on lattice regularization. 

In this context it is interesting to note that the 
one-loop effective potential can be obtained using 

V{H) = Vciassic + const y/V (4) 

where v are the eigenvalues of the second deriva- 
tive of the action with respect to the gauge fields, 
and the summation is over all the eigenvalues. 
This expression is a natural consequence of the 
saddle point approximation: 



V{H) = V,las 



where is the volume factor. The eigenvalue 
of this problem can be found by realizing that 
there is the same symmetry as the one of the Lan- 
dau levels problem, therefore we have exactly the 
same class of solutions. That yield to the eigen- 
values 

^ P + {2n + l + 2S,)gH (6) 

where k = k"^ + kl, n = 0, 1, 2, ... and Sz — ±1 is 
the gluon polarization. 

For the lowest landau mode, i.e. n — 0, Sz — 
— 1, we have v < that give the imaginary part 
of the effective potential. Note that with Abelian 
background it is possible to obtain all eigenvalues 
positive only by insertion of non-gauge invariant 
terms in the action. 



3. LATTICE RESULTS 

Until now we have argued that this lowest Lan- 
dau mode plays a very particular role on the one 
loop perturbative expansion. This motivated us 
to understand what happened to this mode on 
lattice where the task is non-perturbative. 

To extract information on this mode on lattice 
we noted that, due to the finiteness of the lat- 
tice kz, is quantized as well, with kz = 2nnr/Lz, 
where m is an integer. Therefore, the lowest in- 
homogeneous z-mode, m — 1, becomes stable for 
lattices the extent of which in the z-direction is 
smaller than 
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(7) 



This enables us to search for the critical size 
j^criuc^ The homogeneous mode, m = 0, which 
is always unstable, is eliminated by imposing a 
delta condition in the path integral. 

We perform Monte Carlo simulations that gen- 
erate a background field H in the z direction in 
the presence of an Abelian source of strength j. 

We use a heat bath updating procedure with 
periodic boundary conditions and for the com- 
putational technique we address to reference 
To eliminate the homogeneous mode m = we 
force the Polyakov line in the z direction to take a 
fixed value different from zero. Monte Carlo sim- 
ulations are made with a lattice volume * Lz, 
where L is the size of the x,y,t directions. The 
finite effect due to L is not so crucial. In fact, for 
the observable that we are interested, it is suffi- 
cient to use L = 12 — 16. Simulations are made 
by varying Lz from 4 to 50. 

The expectation value of the plaquette in the 
1-2 plane P[Fi2] as a function of (3 and j is mea- 
sured. We monitor the quantity 



X = 



P[Fi2(/3,0)]~P[fi2(/3,j)] 
jP[Fi2(/3,0) 



(8) 



which is proportional to the main contribution 
of the vacuum energy of the plaquette in the z 
direction [^jH , and is very sensible to the presence 
of the unstable mode. The total vacuum energy 
obviously has also contributions from the other 
oriented plaquette; however, in the region we are 
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interested, these other contributions are smoothly 
variable functions, and thus they will not affect 
the critical behavior. 

We evaluated X for different values of (3 and 
j performing, for large j 4500 sweeps after dis- 
charging 500 for thermalization. For smaller j we 
increase the number of sweeps until a significant 
amount of data was collected. 

Our data |2| show that there is no sign of the 
unstable mode away from the critical (3 region 
{P = 2.1 — 2.5). The situation changes dramat- 
ically in the critical region where the instability 
appears as a decrease of the vacuum energy con- 
tribution to the plaquette in the z direction. This 
effect becomes more evident in the presence of 
strong sources. The fact that the instability dis- 
appears outside of the scaling window is a strong 
evidence that the instability is a distinguishing 
feature of the continuum rather than the strong 
coupling vacuum. In order to study the IR prop- 
erty of the continuum theory we should remain 
in the region where the instability is manifest. 
The disappearance of the instability as /? ^ 
may give a clue to understanding the difference 
between the physics of the continuum and the 
strong coupling lattice theory. 

We systematically analyzed the critical region 
of (3 for several values of j. The values of 
L^'^***'^(/3, j) are obtained by interpolating the 
kink of X and taking the median value. 
was found to be dependent on the source j and 
on (3 in this region according to the renormal- 
ization group dependence. It is clear from our 
data that Lf^*'^'^ belongs to the confinement phase 
and that there is good agreement with the renor- 
malization group equation. Our data follow the 
same dependence of the deconfinament transition 
|6| as should be for a dimensional scale length. 
Hence the ratio between the L'^"*'^'^ and the de- 
confinament scale parameter Ldec is independent 
of f3. The relevant physical quantities must be 
obtained in the limit of vanishing of the induced 
field Q{j) 0. From our data we obtain 

j^critic 

= 2.5 ±0.2. (9) 

Using the known value of Ldec @ we can establish 
that the lattice system should be greater than 



Lz ~ 18 in order to reflect the richness of the non- 
perturbative vacuum of the continuum theory. 

In addition to this new scale there is a sequence 
of m modes. In fact for each m we have: 

— ^ = V^ ^ LZ^'^f^—m. (10) 

z{m) "'^^ 

The simulations show that these modes follow 
with good approximation eq. (|lO|) . This is strong 
evidence for the harmonicity of this modes. This 
is a quite surprising result because we are in a re- 
gion of strong non-perturbative effects. It is stim- 
ulating to think in terms of a parallel between the 
situation described above and the integer Quan- 
tum Hall Effect. This similarity is based on the 
dual picture (z <-> i, i? <-> TJ, etc.), and on the 
plateau structures for each mode. 
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